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ABSTRACT 


A  survey  of  the  literature  having  possible  application 
to  the  annular “mist  flow  problem  is  presented. 

The  horizontal  annular-mist  flow  of  natural  gas-water 
mixtures  in  one- and  two  -inch  diameter ,  schedule  40  aluminum 
pipes  was  investigated.  The  input  gas- liquid  volume  ratios 
were  varied  between  10  and  100,000  for  superficial  gas 
velocities  between  25  and  300  ft.  per  sec.  The  mid-point 
pressure  was  held  constant  at  120  psig.  Some  data  for  variable 
mid-point  pressure  and  constant  liquid  and  gas  rates  are  also 
presented  to  show  the  effect  of  gas  density. 

The  pressure  drop  data  for  constant  mid-point  density 
are  presented  as  a  function  of  the  superficial  velocities  of 
the  two  phases.  The  same  data  are  also  presented  in  terms  of 
the  ratio  of  the  two-phase  pressure  gradient  to  the  single - 
phase  gas  pressure  gradient.  Possible  mechanisms  of  the  flow 
in  the  annular-mist  region  are  discussed  and  the  possibility 
of  microscopic  changes  in  flow  pattern  is  indicated.  A 
large  effect  of  a  small  quantity  of  water  was  observed  over 
the  gas  rates  investigated.  Higher  water  rates  had  much  greater 
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effect  at  low  gas  rates  than  at  high  gas  rates.  It  is  pro¬ 
posed  that  the  high  pressure  drop  is  caused  by  the  formation 
of  a  liquid  film  on  the  inside  wall  of  the  pipe.  The  thick¬ 
ness  of  this  film  decreases  as  the  gas  rate  is  increased 
until  it  reaches  some  minimum  value.  The  breakup  of  this 
film  was  not  indicated  over  the  range  of  gas  rates  investigated. 
A  correlation  is  proposed  for  low  liquid  rates  based  on  the 
data  of  this  investigation. 

The 

factor  and 


data  are  also  interpreted  in  terms  of  a  friction 
Reynolds  number  based  on  the  mixture  properties. 
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I  INTRODUCTION 


The  engineer  in  recent  years  has  encountered  many 
problems  concerned  with  the  simultaneous  flow  of  liquids  and 
gases,  immiscible  liquids,  and  solids  and  liquids.  The  petroleum 
engineer  has  been  particularly  concerned  with  the  simultaneous 
flow  of  gases  and  liquids  in  the  reservoir,  well  bore,  gathering 
lines,  transmission  lines  and  processing  plant  equipment.  A 
better  understanding  of  the  two  phase  flow  in  the  well  bore 
could  allow  more  accurate  prediction  of  flowing  bottom  hole 
pressures.  More  recently,  centrally  located  batteries  in  oil 
fields  have  resulted  in  long  gathering  lines  which  must  be 
designed  for  multiphase  flow.  The  design  of  gas-lift  install¬ 
ations  could  also  be  greatly  improved  if  the  flow  pattern  and 
pressure  drop  could  be  predicted  more  accurately. 

A  considerable  amount  of  work,  both  theoretically  and 
experimentally,  has  been  done  towards  a  general  solution  of 
the  problem  of  gas-liquid  flow.  Omer  (1),  Wood  (2)  and  Bennett 
(3)  have  presented  good  literature  surveys  of  the  work  done  on 
two-phase  horizontal  and  vertical  flow. 

The  problem  of  interest  to  the  engineer  involves  the 
prediction  of  the  pressure  drop  and  sometimes  the  flow  pattern 
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from  the  fluid  properties  and  flow  rates  along  with  the  diameter, 
length,  roughness  and  orientation  of  the  pipe. 


Single-phase  Flow 

It  is  instructive  to  first  review  briefly  the  work  that 
has  been  done  on  single-phase  flow.  The  Fanning  friction  factor 
for  flow  in  circular  tubes  is  defined  by  the  equation 

gcp  £Pf 

f  =  2PV2  dL  C1) 


where  f  = 
dPf 

dT~  “ 

V  = 
P  = 
D  - 


Fanning  friction  factor 

pressure  drop  per  unit  length  due  to 
friction  alone 

linear  velocity 

fluid  density 

internal  pipe  diameter 


If  the  extensive  available  data  are  interpreted  in  terms  of 
the  friction  factor  so  defined,  then  for  laminar  flow 


(2) 


where  Re  =  Reynolds  number  for  laminar  flow  (<C  2000). 


This  result  has  also  been  derived  theoretically. 
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For  turbulent  flow  the  value  of  f  cannot  be  predicted 
from  theory  alone,  but  must  be  determined  by  actual  experimental 
work.  By  dimensional  analysis  it  can  be  shown  that,  for 

turbulent  flow,  the  friction  factor  is  a  function  of  the 

DP  V  0 

Reynolds  number,  -  ,  and  of  the  relative  roughness,  ~  , 

r  ^ 

which  is  the  ratio  of  absolute  roughness  to  pipe  diameter. 

Nikuradse  (4),  by  means  of  his  own  experimental  data 
and  that  of  other  investigators  on  smooth  pipes  obtained  the 
following  relationship  between  f  and  Re: 

-pr  =  4.0  log  (Re  y?)  -  0.40  (3) 

vf 

where  Re  >  3000 


Von  Karman  (6),  using  the  universal  velocity-distribution, 
derived  a  theoretical  equation  which  had  the  same  form  as  the 
Nikuradse  equation  and  differed  only  slightly  from  it  in  the 
value  of  the  constants. 


4.06  log  (Re  Jf"  )  -  0.60 


(4) 


The  difference  between  these  two  equations  is  small. 

The  friction  factor  curves  for  rough  pipes  follow  the 
smooth  tube  friction  factor  equation  for  a  short  range  of 
Reynolds  numbers.  Beyond  this  region  there  is  a  transition  to 
fully  developed  turbulent  flow. 
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For  the  transition  region  Colebrook  (7)  proposed  the 
following  empirical  equation: 

V'v 

7=-  “  4  log  2  +  2.28  -  4  log  (1  +  4.67  )  (5) 

where  ^  0.01  (6) 

Re/T 

At  values  of  ^e—  0.01  the  flow  is  fully  turbulent  and 

Re  yj 

the  friction  factor  is  then  independent  of  the  Reynolds  number. 
Von  Karman’s  theoretical  equation  for  fully  developed  turbulent 
flow,  based  on  the  universal  velocity- distribution,  is 

—  =  4  log  —  +  2.28  (7) 

Two -phase  Flow 

In  single -phase  flow  the  energy  loss  resulting  from 
the  irreversibilities  is  due  to  the  viscous  effects  and  eddy 
dissipation.  In  the  simultaneous  flow  of  immiscible  fluids, 
additional  energy  loss  is  caused  by  the  interaction  of  the 
fluids.  The  problem  is  complicated  since  there  are  many  flow 
patterns  and  experience  has  shown  that  the  flow  pattern  has  a 
great  effect  on  the  energy  loss. 

In  a  horizontal  pipe,  when  the  liquid  and  gas  rates  are 
both  small,  the  two  phases  flow  separately  with  liquid  at 
the  bottom  of  the  pipe  and  the  flow  is  called  stratified.  If 
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the  liquid  rate  is  then  held  constant  and  the  gas  rate  is 
increased,  the  liquid  surface  develops  waves  or  ripples.  Further 
increase  in  the  gas  rate  causes  the  liquid  to  climb  the  pipe 
wall  and  to  coalesce  at  the  top.  This  type  of  flow,  which 
consists  of  a  continuous  core  of  gas  surrounded  by  a  layer  of 
liquid  adhering  to  the  tube  wall,  is  called  annular  flow.  More 
rapid  gas  rates  cause  a  portion  of  the  liquid  to  be  dispersed 
into  the  gas  stream  in  the  form  of  mist,  reducing  the  thickness 
of  the  liquid  film.  This  type  of  flow  is  called  annular-mist 
and  is  the  subject  of  the  present  study.  Still  higher  gas 
rates  may  cause  the  liquid  film  to  break  up  and  result  in  a 
mist  flow.  If,  in  the  previously  mentioned  stratified  flow, 
more  liquid  were  flowing  before  the  gas  velocity  was  increased, 
the  liquid  might  fill  the  pipe  for  short  lengths  separated  by 
large  bubbles  of  gas.  In  most  cases  many  small  bubbles  would 
be  mixed  in  the  liquid.  Such  a  flow  is  often  called  slug  flow. 
Still  larger  gas  rates  will  again  cause  annular  and  annular- 
mist  flow  as  in  the  previous  case.  If  the  liquid  and  gas  rates 
were  increased  simultaneously,  the  flow  would  probably  consist 
of  gas  bubbles  dispersed  throughout  the  liquid. 

Any  two  phase  flow  phenomenon  is  defined  uniquely  by  the 
physical  properties  of  each  phase,  the  system  geometry  and  the 
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mass  flow  rate  of  each  phase.  In  reviewing  the  present  state 
of  knowledge  it  is  clear  that  the  relationship  between  the  system 
variables  and  the  energy  loss  due  to  irreversibilities  is  a 
complicated  one.  The  only  hope  of  solving  this  problem  completely 
appears  to  be  in  treating  each  flow  pattern  separately  and  analy¬ 
zing  the  flow  mechanism.  This  point  of  view  is  reflected  in 
the  recent  literature. 

The  present  investigation  is  concerned  with  the  study 
of  horizontal  flow  of  gas -liquid  mixtures  at  high  gas  and 
low  liquid  rates.  In  this  region,  part  of  the  liquid  is  dis¬ 
persed  (or  entrained)  in  the  gas  stream  and  part  of  it  forms 
an  annular  film  on  the  inside  pipe  wall.  This  type  of  mist  or 
annular-mist  flow  is  likely  to  occur  in  condensate  gas  fields 
both  in  the  well  bore  and  the  gathering  lines  when  producing 
at  high  rates. 
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II  LITERATURE  REVIEW 

A  considerable  amount  of  work  has  been  done,  since  the 
1930’s,  on  the  flow  of  gas  liquid  mixtures.  Not  a  great  deal  of 
this  work  applies  to  the  mist  or  annular-mist  flow  regions  and 
the  larger  proportion  of  this  relates  to  vertical  flow.  As 
previously  mentioned  many  investigators  have  presented  good  general 
reviews  of  the  subject  (1,  2,  3).  The  present  review  is  restricted 
to  investigations  with  possible  application  to  the  annular-mist 
flow  region. 

One  of  the  earlier  and  still  very  widely  used  correlations 
is  that  of  Lockhart  and  Martinelli  (7).  They  presented  their 
correlation  as  a  log- log  plot  of  $g  and  $]_  versus  X  with  the 
flow  type  as  parameter,  where 


g 


2 


(8) 


(9) 


and 


X 


2 


(10) 
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The  correlation  was  based  on  data  obtained  in  pipe  sizes 
up  to  one-inch  diameter.  The  values  of  X  for  annular-mist 
flow  are  usually  below  the  range  investigated  by  Martinelli 
(7),  necessitating  extrapolation  of  the  curves.  The  appli¬ 
cation  of  this  correlation  in  the  annular-mist  flow  region 
results  in  calculated  pressure  drops  which  are  too  low  (8) . 

In  addition  the  Martinelli  parameters  have  been  found  inade¬ 
quate  for  the  correlation  of  data  in  the  annular-mist  flow 
region  (8,  9) . 

Chenoweth  and  Martin  (9)  obtained  data  of  low  liquid 

volume  fractions  at  pressures  up  to  100  p.s.i.a.  in  1  1/2- 

inch  and  3-inch  diameter  pipes.  The  deviations  of  the  data 

from  the  Lockhart  and  Martinelli  correlation  were  systematic. 

AP 

A  correlation  has  been  presented  as  a  logarithmic  plot  of  ~ 

versus  liquid  volume  fraction  with  a  parameter  which  reduces 

t°  fg  Pi  for  straight  pipe. 

fi  pg 

Baker  (10*  11)  collected  field  data*  some  of  which 
approaches  the  annular-mist  region*  for  gas-oil  mixtures  in 
4-  to  10-inch  diameter  lines.  The  length  of  lines  varied  from 
2  to  8  miles.  He  found  Martinelli' s  correlation  to  be  inade¬ 
quate  and  presented  equations  for  each  flow  regime  separately. 
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These  equations  employ  Martinelli’s  parameters.  No  equations 


were  presented  for  mist  or  annular-mist  flow, 

Hoogendoorn  (8)  obtained  limited  data  in  the  annular- 
mist  flow  regime  in  24,  50  and  91  mm.  diameter  pipes  with  air 
and  gas-oil  (viscosity  =  2.36  c  .p  . ,  S.G.  =  0.815),  and  also 
in  50  mm.  pipe  with  air  and  spindle  oil  (viscosity  =  20.8 
c.p.,  S.G.  =  0.905).  The  data  for  superficial  gas  Reynolds 
numbers  (Reg)  from  2  X  10^  to  6  X  10^  and  liquid  mass 

velocities  from  9.2  to  855  have  been  present* 

ft?  sec.  _ 

r-.HP. 


have  been  presented  as 

ft?  sec. 

a  plot  of  the  two  phase  friction  factor  f 


(iE)  g 

vdl/  6< 


2  V  pg  - 


based  upon  gas  properties  against  the  superficial  gas  Reynolds 

(Re, 


'g 


D  pg  vg  N 

- - — _  )  number  with  the  water  rate  as  parameter. 


g 


Over  the  limited  range,  Hoogendoorn  has  found  that,  for  liquid 
rates  which  are  not  too  low,  the  friction  factor  is  independent 
of  the  liquid  rate.  Also  no  influence  of  air  density  on 
friction  factor  was  observed  for  air  densities  between  0.075 
and  0.187  lbs.  . 


Govier,  Radford  and  Dunn  (12)  analyzed  the  problem  of 
the  isothermal  flow  of  air  and  water  mixtures  in  vertical  tubes 
and  showed  from  thermodynamic  considerations  that  the  total 
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unit  pressure  drop  could  be  expressed  by  the  relationships: 


VAI/ 


1  +  Rm  +  _i _  2  n  vi 

1  +  1  +  Rv  §c® 


or 


/AP>> 
‘vg  AL} 


_1 

1  + 


1  + 


i  +  si 


R. 


f*  VZ 
g  g 

gcD 


V 


(11) 


(12) 


These  equations  define  a  friction  factor  based  upon 
either  the  liquid  or  gas  properties  as  used  by  Govier  and 
associates.  The  first  term  on  the  right-hand  side  of  the 
equation  is  defined  as  the  hydrostatic  head  component  and  the 
second  term  as  the  irreversibility  component.  Actually, 
however,  the  first  term  represents  the  hydrostatic  head  only 
if  there  is  no  holdup.  Govier  included  the  entire  effect  of 
holdup  in  the  second  term.  It  is,  therefore,  not  possible 
to  compare  the  irreversibility  component  from  this  analysis 
with  the  energy  loss  due  to  irreversibility  in  horizontal 
flow. 

Bennett  and  Thornton  (13)  calculated  the  friction  fac- 
tor  based  on  liquid  properties  as  defined  by  Govier  et  al. 

The  experiments  were  conducted  in  vertical  glass  pipes  of 
1.36 -inch  and  1.25- inch  internal  diameter  and  also  in  an 
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annulus  between  1-inch  outside  diameter  and  1.5- inch  inside 
diameter  glass  pipes.  This  annular-mist  flow  data  has  been 
graphically  represented  by  two  straight  lines  on  a  log-log 
plot  where  f-j  x  ReCT  is  plotted  against  Re-i .  The  data  for  the 
annulus  and  the  larger  pipe  are  represented  by  one  line  and 
the  data  for  the  smaller  pipe  by  the  other  line.  Slopes  of 
these  lines  are  -2.81  and  -2.75  respectively. 

In  a  recent  study  Ros  (14)  has  made  use  of  dimensional 
analysis  to  study  the  simultaneous  flow  of  oil  and  gas  in 
vertical  pipes.  He  has  derived  an  equation  for"mist  flow"  which 
is  only  valid  for  flow  at  atmospheric  conditions.  This  type 
of  study  shows  promise  in  multiphase  flow  applications. 

Teletov  (15)  based  his  correlation  on  a  friction  factor 
and  Reynolds  number  based  on  the  mixture  properties  which  may 
be  represented  as 


Sc  » 


tpm  2(Vl  +  Vg)2  Pm 


Re. 


m 


D  Pi  V! 

M-l 


dP 

dL 


D  Pg  Vg 


M- 


g 


(13) 

(14) 


The  correlation  was  presented  in  terms  of  a  deviation 
factor*  Y  *  as  a  function  of  the  volume  fraction  of  liquid 
and  the  viscosity  ratio  of  the  two  phases.  The  deviation  factor 
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is  defined  as 


where 


Y 


f 

•Ltpm 

f* 


(15) 


f*  =  friction  factor  at  Rem  from  single- 
phase  charts  or  equations. 

The  object  of  the  correlations  of  Lockhart  and  Mar- 
tinelli  (7)*  Govier  et  al  (12),  Bennett  and  Thornton  (13), 
Chenoweth  and  Martin  (9)  and  others  of  this  type  is  to  allow 
the  prediction  of  pressure  drop  directly  from  the  flow  vari¬ 
ables  without  requiring  a  mechanistic  analysis  of  the  flows 
in  the  various  regions.  The  total  energy  dissipation  due  to 
irreversibilities  is  usually  expressed  in  terms  of  a  two-phase 
friction  factor  or  some  other  single  factor. 

For  the  specific  case  of  annular°mist  flow,  where 
measurements  of  entrainment  and  film  thickness  are  available, 
a  different  type  of  analysis  may  be  applied.  In  this  type 
of  analysis  the  flow  in  the  liquid  film  and  that  in  the  gas 
core  with  entrained  liquid  are  treated  separately  to  calculate 
the  pressure  drop.  This  general  approach  has  been  taken  by 
Wicks  and  Dukler  (16),  Calvert  and  Williams  (17),  Anderson  and 
Mantzouranis  (19)  and  Collier  and  Hewitt  (20)  and  McManus  (21). 
An  effort  has  been  made  by  these  investigators  to  understand 
the  flow  mechanism  by  a  study  of  the  microscopic  nature  of  the 
flow  pattern. 
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Wicks  and  Dukler  (16)  have  reported  entrainment  data 
for  annular-mist  flow  of  air-water  mixtures  in  one-and  three- 
inch  diameter  horizontal  pipes.  The  pressure  drop  data  were 
presented  only  for  the  one-inch  pipe.  Their  investigation 
covered  air  rates  from  125  to  475  lbs. /hr.  and  water  rates 
from  521  to  3126  lbs. /hr.  The  influence  of  entrance  section 
design  on  pressure  drop  and  entrainment  was  demonstrated.  A 
correlation  relating  an  " entrainment  variable”  to  the  ratio 
of  the  single-phase  pressure  drop  gradients  was  also  presented. 

The  analyses  of  Calvert  and  Mantzouranis  were  derived 
for  upwards  vertical  flow.  These  have  been  applied  to  the 
flow  data  of  Bennett  and  Thornton  and  others  by  Collier. 

The  object  of  the  analysis  of  Calvert  and  Williams  (17) 
was  to  calculate  the  liquid  flow  in  an  annular  film  from  a 
knowledge  of  the  film  thickness  and  the  pressure  drop.  This 
was  achieved  by  postulating  a  velocity  distribution  for  the 
film  in  terms  of  shear  (and  thus  pressure  drop) .  The  film  was 
assumed  to  consist  of  a  laminar  layer  adjacent  to  the  wall 
and  a  turbulent  layer.  The  velocity  distributions  for  these 
two  layers  were  integrated  with  respect  to  the  distance  from 
the  wall  to  yield  the  total  liquid  flow  rate. 

Anderson  and  Mantzouranis  (18,  19)  have  used  the 
universal  velocity  profile  of  von  Karman.  Integration  of  this 
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velocity  profile  to  yield  film  flow  rates  from  film  thickness 
and  pressure  drop  data  might  be  more  accurate  than  the  Calvert 
derivation  since  the  former  takes  account  of  the  existence  of 
the  "buffer"  region  between  the  laminar  and  turbulent  regions. 

Collier  and  Hewitt  (20)  applied  the  Mantzouranis  and 
the  Calvert  analyses  to  predict  the  liquid  film  flow  rates. 
Mantzouranis'  analysis  predicted  the  rates  more  closely  than 
the  Calvert  analysis.  At  high  gas  rates  and  low  liquid  rates 
the  predicted  film  flow  rates  were  much  higher  than  the  actual 
in  both  the  Calvert  and  Mantzouranis  analyses.  Two  possible 
reasons  for  deviation  are  givens 

(1)  Interchange  of  droplets  between  the  gas  core 

and  the  liquid  film  -  the  resulting  acceleration/ 
deceleration  effects  causing  a  loss  of  energy. 

(2)  The  collapse  of  surface  waves  on  the  liquid 
film  with  consequent  dissipation  of  energy. 

Since  the  wave  system  is  maintained  this  would 
result  in  pressure  loss.  The  waves  might  also 
have  a  profound  effect  on  the  velocity  distribu¬ 
tion  in  the  film. 

McManus  (19)  has  reported  on  the  local  values  of  film 
thickness  and  wave  heights  in  horizontal  two-phase  flow.  He 
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has  shown  that  appreciable  eccentricity  of  the  film  exists 
and  that  the  wave  height  is  exponentially  related  to  the  film 
depth.  It  has  been  suggested  that  these  effects  must  be 
accounted  for  in  any  complete  analysis  of  horizontal  flow. 
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III  EXPERIMENTAL  EQUIPMENT 

The  equipment  was  designed  to  obtain  pressure  drop 
data  for  a  natural  gas-water  system  at  high  superficial  gas 
velocities  (30  to  300  ft/sec.)  and  low  superficial  water 
velocities  (0,001  to  1  ft/sec.).  The  flow  control  and  measure¬ 
ment  equipment  was  selected  to  give  the  maximum  desired  velocities 
in  a  one-inch  diameter  pipe.  The  range  of  velocities  investigated 
for  the  two -inch  pipe  was  lower  with  the  same  maximum  mass  flow 
rates  as  in  the  one-inch  pipe. 

Arrangements  were  made  with  Northwestern  Utilities  Limited 
(the  natural  gas  distribution  company  for  the  city  of  Edmonton) 
to  install  all  equipment  at  their  City  Gate  Station  No.  1. 

Natural  gas  at  sufficiently  high  pressure  was  available  at 
this  place. 

Two  test  sections  of  aluminum  schedule  40  pipe  having 
internal  diameters  of  1.049  and  2.067  inches  were  installed 
as  shown  in  Figure  1.  The  one-inch  test  section  was  40  feet 
long  with  approximately  15  feet  of  calming  section  upstream  and 
downstream  of  the  test  section.  Five  pressure  taps  spaced  10 
feet  apart  were  provided  on  the  one-inch  test  section.  The  two- 
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inch  test  section  was  80  feet  long  with  approximately  15  feet 
of  the  calming  section  upstream  and  downstream  of  the  test 
section.  Five  pressure  taps  spaced  20  feet  apart  were  provided 
on  the  two-inch  test  section.  The  pipe  was  purchased  in  20-foot 
long  sections.  A  welded  joint  is  shown  in  Figure  2.  Detail 
of  the  pressure  taps  is  shown  in  Figure  3.  Particular  care 
was  taken  to  keep  the  pipe  smooth  both  at  the  joints  and  at 
the  pressure  taps. 

Gas  System 

Natural  gas  from  the  gas  plant  at  Pembina  oil  field  in 
Alberta  (average  composition  given  in  Appendix  A-l)  was  available 
at  the  City  Gate  Station  at  a  variable  pressure  of  300  to  600 
p.s.i.g.  The  gas  was  passed  through  a  pressure  control  valve, 
two  electric  heat  exchangers  in  series  and  an  orifice  flow 
meter  to  a  tee  where  it  could  be  directed  to  either  the  one-inch 
or  the  two-inch  test  section.  Water  was  injected  axially  into 
the  gas  stream  at  the  beginning  of  the  calming  section,  for 
each  section.  Spray  nozzles  were  used  to  spray  the  water  into 
the  gas  stream.  Gas  and  water  from  the  end  of  the  test  section 
were  taken  to  a  flare  stack  approximately  300  feet  from  the  end 
of  the  test  section.  The  flare  stack  was  made  of  a  40  foot 
vertical  section  of  four-inch  pipe  which  was  located  about  200 
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feet  from  the  City  Gate  Station  facilities .  At  the  flare 
stack  gas  was  burned  and  water  was  either  evaporated  or  sprayed 
around  the  flare  as  the  gas  burned.  During  the  dry  gas 
calibration  runs  the  gas  was  returned  to  a  lower  pressure  gas 
line  operating  at  50  p.s.i.g. 

The  gas  flow  rate  was  regulated  by  a  two-inch,  single¬ 
port  pneumatic  valve  with  3/4-inch  inner  valve  (Fisher  667-D). 

The  control  valve  was  operated  by  a  3  to  15  p.s.i.g.  air  signal 
to  the  bottom  of  the  diaphragm  from  a  proportional  plus  reset 
controller.  The  input  to  the  controller  was  from  a  pressure- 
transmitter  installed  at  the  mid-point  of  either  one  of  the 
two  test  sections.  Mid-point  pressure  was  recorded  on  a  chart 
to  observe  the  fluctuations,  if  any,  in  the  pressure.  The 
temperature  of  the  gas  dropped  to  as  low  as  10°F  during  expansion 
through  the  control  valve.  It  was  then  heated  to  above  the 
freezing  point  of  water  and  as  close  to  the  ambient  temperature 
as  possible  by  two  electric  heat  exchangers  (Chromolox  Type 
GCH-6120W) .  The  heat  exchangers  were  of  3  K.W.  and  8  K.W. 
capacity  with  on-off  type  controllers. 

The  heated  gas  was  metered  by  a  two-inch  orifice  meter 
constructed  to  American  Gas  Association  specifications.  The  orifices 
fitting  used  (Senior  Daniel  Orifice  Fitting)  allowed  the  changing 
of  the  orifice  under  pressure.  The  static  and  differential 
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pressures  were  recorded  on  a  mercury  type  flow  recorder.  The  gas 
temperature  was  measured  downstream  of  the  orifice  by  a  dial 
thermometer . 

Water  System 

City  water  from  300- gallon  water  tanks  was  pumped  by 
two  centrifugal  pumps*  which  could  be  used  either  individually 
or  in  series.  The  pumps  had  3  H.P.  and  5  H.P.  motors  respectively. 
The  water  rate  was  measured  and  controlled  by  an  integral 
orifice  differential  pressure  (D/P)  cell  (Foxboro  Type  13  A)* 
proportional  plus  reset  rate  controller  recorder  and  two  control 
valves.  The  control  valves  of  1/8-inch  and  1/2-inch  orifice 
size  were  installed  in  parallel  upstream  of  the  integral  ori¬ 
fice  D/P  cell  so  that  each  valve  could  be  used  independently. 

Two  valves  were  required  because  neither  individually  could 
give  adequate  control  over  the  range  of  water  rates  desired. 

Two  lines  with  independent  valves  were  provided  from  the  down¬ 
stream  of  the  D/P  cell  to  the  upstream  of  each  one  of  the  two 
test  sections.  The  arrangement  of  the  equipment  is  shown  in 
Figure  1. 
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Pressure  Drop  Measurement 

The  pressure  drop  in  the  one-inch  pipe  was  measured 
by  an  electrical  D/P  cell  (Foxboro  Resistance  Dynalog  2400 
Series)  having  a  range  of  0-600  "  H2O  pressure  drop  and  a 
multipoint  recorder  (Foxboro  Dynalog) .  The  pressure  taps  are 
numbered  1  to  5  starting  at  the  upstream  side.  The  manifold 
arrangement  for  the  D/P  cell  is  shown  in  Figure  4.  At  high 
flow  rates  the  pressure  drops  between  pressure  points  1-2* 

2-3*  3-4  and  4-5  were  measured  while  at  low  flow  rates  the 
pressure  drops  between  points  1-2*  1-3*  1-4  and  1-5  were  meas¬ 
ured.  This  allowed  the  measurement  of  pressure  drops  over 
the  test  section  of  up  to  70  p.s.i.  with  a  maximum  estimated 
error  of  ±  6"  H2O.  The  pressure  drop  in  the  two-inch  pipe  was 
measured  by  two  pneumatic  D/P  cells  having  range  of  0-100  inch 
H2O  pressure  drop  (Minneapolis  Honeywell  Brown  Converters)  and 
a  recorder  (Taylor  Fulscope) .  The  arrangement  is  shown  in 
Figure  5.  The  pressure  drops  1-3  and  3-5  were  measured.  At 
very  high  flow  rates  only  the  pressure  drop  between  points 
2-3  was  measured. 

A  transparent  section  of  lucite  pipe  6M  long  was  inserted 
downstream  of  the  one-inch  test  section  to  observe  and  photo¬ 
graph  the  flow  pattern.  Plate  I  shows  the  transparent  section 
in  position. 
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PLATE  X 


TRANSPARENT  SECTION  IN  ONE-INCH  PIPE 
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Instrument  Air  Supply 

For  the  operation  of  all  pneumatic  instruments  air  at  a 
pressure  of  20  p.s.i.g.  was  required.  This  air  was  supplied 
by  an  air-compressor  through  a  surge  tank  at  100  p.s.i.g. 

The  pressure  was  reduced  to  40  p.s.i.g.  at  the  outlet  of  the 
surge  tank  and  then  to  20  p.s.i.g.  near  the  instruments. 

The  two- stage  reduction  of  pressure  and  the  surge  tank  elimina¬ 
ted  the  pulsations  caused  by  the  reciprocating  compressor. 

Plate  II  shows  the  compressor,  air  filter,  surge  tank  and  the 
first- stage  pressure  regulator. 

Power  Supply 

Electrical  power  was  supplied  to  the  motors  and  the 
heat  exchangers  from  a  3-phase,  220-volt  switchboard.  Overload 
protection  switches  were  provided  for  each  motor  separately. 
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PLATE  II 

INSTRUMENT  AIR  SUPPLY 
AIR  COMPRESSOR  AND  SURGE  TANK 
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IV  EXPERIMENTAL  PROCEDURE 

The  pressure  drop  and  flow  rate  measuring  equipment 
were  first  calibrated  as  outlined  in  Appendix  C.  The  single¬ 
phase  pressure  drop  for  gas  Reynolds  numbers  of  approximately 
3  x  10^  to  2.5  x  10^  was  then  obtained  to  determine  the 
relative  roughness  of  the  test  sections.  The  single-phase 
data  are  presented  in  Figures  C-2  and  C-3  of  Appendix  C. 

Experiments  were  then  conducted  to  obtain  pressure 
drop  and  flow  rate  data  for  the  simultaneous  flow  of  gas  and 
water  at  a  constant  mid-point  pressure  and  temperature.  A 
mid-point  pressure  of  120  p.s.i.g.  was  selected.  Actually  it 
was  only  possible  to  maintain  the  mid-point  pressure  between 
118  and  121  p.s.i.g.  and  the  average  temperature  between  36°F 
and  73°F.  Most  of  the  data,  however,  was  obtained  at  a  constant 
mid-point  pressure  of  120  p.s.i.g.  and  an  average  temperature 
of  about  50°F.  The  superficial  water  velocities  were  varied 
for  the  one-inch  pipe  from  0.0067  ft. /sec.  to  2.42  ft. /sec. 
and  for  the  two-inch  pipe  from  0.00303  ft. /sec.  to  0.525  ft. /sec. 
The  superficial  gas  velocities  were  varied  for  the  one-inch 
pipe  from  40  ft. /sec.  to  300  ft. /sec.  and  for  the  two-inch  pipe 
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from  25  ft. /sec.  to  130  ft. /sec.  The  pressure  drop  data  were 
obtained  at  constant  superficial  water  velocities  by  varying 
the  superficial  gas  velocities  over  the  range.  Water  required 
to  saturate  the  gas  for  any  gas  velocity  was  calculated  and 
added  to  the  water  required  to  maintain  a  certain  constant 
superficial  water  velocity  in  the  pipe  line.  This  water  rate 
was  then  set  on  the  water  flow  controller.  For  high  water 
rates  the  amount  of  water  required  to  saturate  the  gas  was 
negligible.  The  variables  measured  are  presented  in  Table  1. 
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TABLE  1 
DATA  OBTAINED 

GAS  FLOW  MEASUREMENT  (2M  Orifice  Meter) 

Orifice  size  in  inches 

Static  pressure  (Downstream)  in  p.s.i.g. 
Differential  pressure  in  inches  of  water 
Gas  temperature  in  °F 
Atmospheric  pressure  in  p.s.i.a. 

WATER  FLOW  MEASUREMENT  (Integral  Orifice  D/P  cell) 
Orifice  size  in  inches 

Differential  pressure  in  inches  of  water 

ONE -INCH  TEST  PIPE  (1.049-inch  internal  diameter) 

Pressure  drops  between  pressure  taps  1-2,  2-3, 

3-4  and  4-5  in  inches  of  water 

Mid-point  pressure  (at  pressure  tap  3)  in  p.s.i.g. 

Upstream  temperature  in  °F 

Downstream  temperature  in  °F 

A 

TWO- INCH  TEST  PIPE  (1.067- inch  internal  diameter) 

Pressure  drop  between  pressure  taps  1-3  and  3-5  in 
inches  of  water 
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Two- Inch  Test  Pipe  (Continued) 

Mid-point  pressure  (at  pressure  tap  3  )  in  p.s.i.g. 
Upstream  temperature  in  °F 
Downstream  temperature  in  °F 


WATER  CONTENT  OF  THE  GAS 


Weight  of  the  desiccant  tube  before  gas  flow  in  grams 

Gas  flow  through  the  desiccant  in  cubic  feet  at  the 
temperature  of  gas  and  2  to  6  inches  of  water 
above  atmospheric  pressure 

Weight  of  the  desiccant  tube  after  the  measured  gas 
flow  through  the  desiccant  in  grams 
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Experimental  data  at  approximately  constant  superficial 
gas  velocity,  constant  superficial  water  velocity  and  variable 
mid-point  pressure  was  obtained  by  establishing  critical  flow 
conditions  at  the  downstream  back  pressure  control  valve.  When 
the  pressure  upstream  of  the  valve  was  about  two  times  or  more 
the  pressure  downstream  of  the  valve,  the  superficial  gas 
velocity  was  constant. 
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V  EXPERIMENTAL  RESULTS 

The  composition  and  properties  of  the  gas  and  the 
geometry  of  the  system  are  presented  in  Tables  A-l  and  A-2 
of  Appendix  A.  Appendix  B  contains  all  of  the  pressure  drop 
data.  Section  B-9  of  Appendix  B  describes  the  calculation 
procedure  for  the  calculated  data  in  this  Appendix.  The 
calibration  data  and  other  information  regarding  the  cali¬ 
bration  of  the  equipment  are  included  in  Appendix  C. 


Single-phase 

The  results  of  the  tests  with  gas  flowing  alone  are 
given  in  Appendix  B  in  Tables  B-l  and  B-3  for  the  one-inch 
pipe  and  Tables  B-2  and  B-4  for  the  two-inch  pipe.  The  inte¬ 
grated  form  of  the  Fanning  equation  (16)  for  gas  flow  was 
used  to  calculate  the  friction  factor  for  each  of  the  four 
10- foot  long  sections  of  the  one-inch  test  pipe.  Since  for 
certain  of  the  tests  the  Reynolds  numbers  were  sufficiently 
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where  Q  = 


L  = 
d  = 


f  = 


gas  flow  measured  at  TQ  and  P  ,  std.  cu. 

ft.  per  hr. 

length  of  line,  miles 

internal  diameter  (I.D.),  in. 

upstream  pressure 

downstream  pressure 

gas  gravity  (air  =  1.0) 

average  compressibility  factor 

Fanning  friction  factor 


high  that  fully  developed  turbulence  was  to  be  expected.  The 
friction  factor  for  these  tests  was  assumed  independent  of 
Reynolds  number  and  is  a  function  only  of  the  relative  rough¬ 
ness.  The  equation  of  Nikuradse  was  then  used  to  calculate 
the  relative  roughness  for  each  section.  The  single-phase 
data  are  presented  for  each  of  the  sections  of  the  one-inch 
pipe  in  Figure  C-2  of  Appendix  C.  The  average  relative  rough¬ 
nesses  for  the  various  sections  of  the  one-inch  pipe  were 
0.00006,  0.00009,  0.00015  and  0.0002. 

An  average  value  for  the  mid-point  relative  roughness 
was  then  obtained  by  calculating  the  slope  of  the  pressure 
versus  length  curves  at  this  point  for  all  of  the  flow  rates. 
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This  slope  was  then  used  to  calculate  the  friction  factors 
by  the  differential  form  of  the  Fanning  equation.  A  similar 
procedure  was  followed  for  the  two-inch  pipe.  The  friction 
factor  versus  Reynolds  number  curves  for  the  two  pipes  are 
presented  in  Figure  C-3  of  Appendix  C.  The  data  for  the  one- 
inch  pipe  are  best  represented  by  the  Colebrook  equation  for 

«  0.00009  and  the  data  for  the  two-inch  pipe  by  the  equation 

for  —  =  0.000048.  The  Colebrook  curves  for  these  two 

D 

relative  roughnesses  are  also  shown  on  Figure  C-3.  The  rela¬ 
tive  roughness  for  the  two -inch  pipe  is  approximately  one 

half  that  of  the  one-inch  pipe.  This  gives  approximately 

; 

I 

the  same  absolute  roughness  of  0.00009  inches  for  the  two 
sizes  of  aluminum  pipe.  This  is  a  reasonable  value  for  the 
roughness  of  aluminum  pipe*  confirming  the  technique.  A 
good  check  was  also  obtained  between  the  calculated  and  the 
experimental  data  for  single-phase  flow  as  shown  in  Figures  7 
and  8. 

Two-phase  Data 

The  experimental  and  calculated  two-phase  data  are 
tabulated  in  Appendix  B.  Tables  B-l  and  B-2  contain  the 
experimental  data  for  the  one- and  two-inch  pipes  respectively. 
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while  tables  B-3  and  B-4  contain  the  calculated  data.  The 
calculated  data  in  these  tables  consist  of  the  gas-phase  density 
at  the  mid-point,  the  superficial  velocities  of  the  two  phases, 
the  mid-point  pressure  gradient  and  the  mid-point  pressure 
gradient  adjusted  to  a  constant  gas  density  of  0.47  lbs.  per 


The  two-phase  pressure  gradient  at  the  mid-point  of 
dP 

the  test  section,  (■— -)  ,  was  also  obtained  by  measuring 

dL 

the  slope  of  the  pressure  versus  length  curves  at  this  point 
for  each  test.  Typical  pressure  gradient  curves  for  the  one- 
inch  pipe  are  shown  in  Figure  6  for  various  flow  rates. 

Figures  7  and  8  and  Tables  B-3  and  B-4  represent  the 
basic  pressure  drop  data  adjusted  to  a  constant  gas  density 
for  the  two  pipes.  The  figures  show  the  pressure  gradient 
plotted  against  the  superficial  gas  velocity  for  various  con¬ 
stant  superficial  water  velocities.  The  single-phase  calcula¬ 
ted  lines  and  experimental  data  points  are  also  shown  in 
Figures  7  and  8  for  comparative  purposes.  The  single-phase 
lines  were  calculated  by  assuming  relative  roughnesses  of 
0.00009  and  0.000048  for  the  one- and  two-inch  pipes  respectively. 

Cross  plots  of  the  data  of  Figures  7  and  8  are  presented 
in  Figures  9  and  10  respectively.  These  show  the  pressure 
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FIGURE  6 


TYPICAL  ONE-INCH  PRESSURE  GRADIENT 
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gradient  plotted  against  the  superficial  water  velocity  for 
various  constant  superficial  gas  velocities.  The  data  are 
extrapolated  to  the  appropriate  single-phase  gas  pressure 
gradient  values  (the  addition  of  the  small  constant  to  the 
superficial  water  velocity  permitted  the  representation  of 
zero  water  velocity  on  logarithmic  coordinates) .  Lines  of 
constant  gas -liquid  volume  ratio  are  also  shown  on  these 
figures . 

From  Figures  9  and  IQ  smoothed  pressure  drop  data  at 

selected  superficial  water  velocities  were  obtained.  From 

these  smoothed  data  the  ratio  of  the  two-phase  pressure 

dP  HP 

gradient  to  the  single-phase  pressure  gradient,  (—)/(—)  , 

dL  dL  g 

was  calculated.  The  calculated  data  are  tabulated  in  Table 
B-5  and  B-6  of  Appendix  B  for  the  one- and  two-inch  diameter 
pipes  respectively.  Superficial  water  velocities  of  0.01, 
0.03,  0.10,  0.30,  1.0  and  3.0  ft.  per  sec.  were  selected  for 
the  one-inch  pipe  and  values  of  0.001,  0.003,  0.010,  0.03, 
0.10,  0.30  and  1.0  ft.  per  sec.  were  taken  for  the  two-inch 
pipe.  The  highest  velocity  in  each  case  is  an  extrapolation 
of  the  experimental  data.  Superficial  gas  velocities  for  the 
one- inch  pipe  ranged  from  30  to  300  ft.  per  sec.  and  those  for 
the  two-inch  pipe  from  25  to  120  ft.  per  sec. 
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The  data  of  Tables  B-5  and  B-6  are  represented  by 

Fgiures  11  and  12  where  (4r)/(4r)  is  plotted  against  the 

dL  dL  g 

superficial  gas  velocity  for  various  superficial  water 
velocities. 

The  data  were  also  obtained  from  Figures  9  and  10 
at  constant  gas-liquid  volume  ratios  and  are  tabulated  in 
Tables  B-7  and  B-8  of  the  Appendix  B.  For  the  one-inch  pipe, 
gas  to  liquid  volume  ratios  of  50,  100,  200,  400,  1000,  2000, 
4000,  10,000,  and  100,000  and  for  the  two-inch  pipe  40,  100, 
200,  400,  1000,  2000,  4000,  10,000,  20,000,  40,000,  100,000 
and  1,000,000  were  selected.  Figures  13  and  14  represent  the 
data  of  Tables  B-7  and  B-8  where  the  friction  factor  based 
on  mixture  properties  is  plotted  against  the  Reynolds  number 
of  the  mixture  for  constant  gas-liquid  volume  ratios. 

Limited  data  at  variable  mid-point  density  and  constant 
superficial  velocities  of  the  two-phases  were  also  obtained. 
These  data  are  included  in  Tables  B-l,  B-2,  B-3  and  B-4  of 
Appendix  B  along  with  the  constant  density  data.  Figure  15 
indicates  the  variation  of  the  pressure  drop  ratio  (two-phase 
to  single-phase)  with  the  gas  density  at  approximately  constant 
flow  rates  of  the  phases. 
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D  40  60  80  100  200  300 

SUPERFICIAL  GAS  VELOCITY,  Vg  ,  ft.  per  sec. 


TWO-INCH  DIAMETER  PIPE. 
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SUPERFICIAL  GAS  VELOCITY 


CONSTANT  GAS  TO  WATER 
- VOLUME  RATIO,  Rv  =  200 
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4  X  I05  6  X  I05  10s 

REYNOLDS  NUMBER  Rem,  BASED  ON  MIXTURE  PROPERTIES. 
FIGURE  13,  ONE- INCH  PIPE. 
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Observation  of  Flow  Pattern 

The  flow  pattern  was  observed  through  a  transparent 
section  of  one-inch  pipe  over  the  range  of  flow  rates  inves¬ 
tigated.  Visual  observations  showed  a  milky  colour  mist  over 
the  entire  range  of  flow  rates.  At  high  water  rates  and  low 
gas  rates  there  was  some  indication  of  surface  waves 
on  the  inside  pipe  wall.  Considerable  difficulty  was  experienced 
in  photographing  the  flow  pattern.  A  35  mm.  camera  with 
shutter  speeds  as  fast  as  1/1000  sec.  and  a  movie  camera  with 
speeds  of  up  to  8000  frames  per  sec.  (effective  exposure  per 
frame  of  1/24*000  sec.)  were  used.  Three  photographs  of  the 
flow  pattern  are  included  in  Plate  III.  Photograph  1  taken 
at  a  high  water  rate  and  low  gas  rate  shows  the  surface  waves 
in  the  liquid  film  on  the  inside  pipe  wall.  At  approximately 
the  same  gas  velocity  and  a  low  liquid  rate  Ihotograph  2 
shows  no  surface  waves.  At  high  gas  rates  and  over  the  entire 
range  of  water  rates  photographic  observations  of  the  flow 
pattern  were  similar  to  Photograph  3.  There  is  no  indication 
of  surface  waves  in  Photograph  3.  This*  however*  does  not 
mean  a  liquid  film  was  not  present.  The  location  of  the  test 
equipment  made  it  difficult  to  properly  light  the  flow  for 
photographic  purposes.  Better  photographic  techniques  might 
show  the  liquid  film  that  probably  is  present  at  these  velocities. 
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Photograph  1 


Va  =  60  ft. /sec. j  Vw 
§ 


2.42  ft. /sec. 


Photograph  2  Vg  =  60  ft. /sec. ,  Vw  =  0.4  ft. /sec. 


Photograph  3  V  =  300  ft. /sec.,  Vw  =  0,4  ft. /sec. 
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PLATE  III 
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VI  DISCUSSION  OF  RESULTS 


The  experiments  have  indicated  a  large  effect  of  a 

.  ....  / 

small  quantity  of  water  on  the  pressure  drop  in  the  pipeline 
flow  of  gas  at  high  velocities.  This  effect  may  best  be 
shown  by  plotting  the  ratio  of  the  two-phase  pressure 
gradient  to  the  single-phase  pressure  gradient  as  a  function 
of  the  superficial  gas  velocity.  In  this  type  of  represen¬ 
tation  either  gas-liquid  volume  ratio  or  the  liquid  rate  may 
be  the  parameter.  Figures  11  and  12  are  representations  of 
this  type  for  the  one- and  two-inch  diameter  pipes  with  the 
water  rate  as  the  parameter.  The  diameter  and  the  roughness 
of  the  pipe  have  different  effects  on  the  flow  pattern,  and 
hence  on  the  pressure  drop,  at  different  flow  rates.  The 
pressure  gradient  ratio  for  high  water  rates  drops  sharply  as 
the  gas  rate  is  increased.  This  same  effect,  however,  is  not 
observed  for  the  two-inch  pipe.  At  a  superficial  gas  velocity 

of  30  ft.  per  sec.  and  a  superficial  water  velocity  of  1.0 

dP  dP 

ft.  per  sec.j  the  pressure  gradient  ratio,  —  /(■—),  f°r 

the  one-inch  diameter  pipe  is  9.5  and  for  the  two-inch  pipe, 
4.8.  At  a  higher  superficial  gas  velocity  of  100  ft.  per  sec. 
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and  the  same  water  rate  the  pressure  gradient  ratio  for  the 
one- inch  pipe  drops  to  4.7  and  for  the  two -inch  pipe  to  4.5. 

The  pressure  gradient  ratio  curves  for  the  two-inch 
diameter  pipe  exhibit  maximum  and  minimum  values  and  behave 
similarly  over  the  range  of  flow  rates  investigated.  At  lower 
water  flow  rates  the  effect  of  diameter  on  the  flow  pattern 
seems  to  be  small  and  the  pressure  gradient  ratio  curves 
for  the  one-inch  pipe  are  similar  in  shape  to  the  curves  for 
the  two-inch  pipe.  Below  a  superficial  gas  velocity  of  approxi¬ 
mately  200  ft.  per  sec.  the  pressure  gradient  ratio  in  the 
one-inch  diameter  pipe  decreases  with  increasing  gas  rate. 

Further  increase  in  gas  velocity  to  300  ft.  per  sec.  causes 
little  change  in  the  pressure  gradient  ratio,  especially  for 
low  water  rates.  The  maximum  superficial  gas  velocity 
investigated  in  the  two-inch  diameter  pipe  was  120  ft.  per 
sec.  Figure  12  for  the  two-inch  pipe  shows  a  fairly  sharp 
drop  in  the  pressure  gradient  ratio  for  high  water  rates  as 
the  gas  velocity  is  increased  from  about  60  to  120  ft.  per 
sec.  The  data  indicate  that  at  higher  velocities  the  effect 
of  water  rate  on  pressure  drop  decreases  in  the  two-inch  diameter 
pipe  in  a  manner  similar  to  the  one- inch  diameter  pipe. 

Although  all  of  the  data  are  in  the  annular-mist  flow 
region,  the  changes  in  shape  of  the  curves  may  be  interpreted 
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as  representing  the  microscopic  changes  in  flow  pattern. 

This  is  particularly  clear  from  the  two-inch  data  as  shown  in 
Figure  12.  The  loci  of  the  maxima  and  the  minima  of  these 
curves  divide  the  data  into  three  regions.  The  changes  in 
film  thickness  and  the  film  surface  characteristics  are  probably 
the  more  important  factors  influencing  the  shape  of  these 
curves. 

The  single-phase  pressure  drop  lines,  as  shown  by 
Figures  7  and  9  are  almost  straight  and  parallel  for  the  two 
test  pipes  over  the  range  of  gas  velocities  investigated.  The 
difference  in  the  diameter  and  roughness  effect  between  the 
two  pipes  may  be  eliminated  for  single-phase  flow  by  multiplying 
the  two-inch  pipe  pressure  gradient  by  2.30.  The  factor  2.30 
is  the  ratio  of  the  pressure  gradients  for  the  two  pipes  in 
single-phase  gas  flow.  Figure  16  shows  the  actual  pressure 
gradient  for  the  one-inch  pipe  and  the  pressure  gradient 
multiplied  by  2.30  for  the  two-inch  pipe  at  superficial  water 
velocities  of  0,  0.01,  0.1  and  1.0  ft.  per  sec.  for  both 
pipes.  The  agreement  between  the  data  of  two  pipe  sizes  is 
good  for  water  velocities  of  0.01  and  0.1  ft.  per  sec.  The 
super imposition  of  the  curves  at  low  water  rates  shows  the 
small  effect  of  diameter  on  flow  pattern.  At  a  higher  water 
velocity  of  1.0  ft.  per  sec.  the  effect  of  diameter  increases 
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and  the  data  cannot  be  superimposed  in  the  same  manner.  The 
data  for  the  two-inch  pipe  at  a  superficial  water  velocity 
of  1.0  ft.  per  sec.  are  extrapolated  from  Figure  10.  This 
correlation  may  be  used  for  estimating  two-phase  pressure 
drops  in  the  annular-mist  flow  region  in  small  diameter  pipes 
at  low  water  flow  rates. 

For  a  generalized  solution  of  the  two-phase  flow  problem 
it  is  advantageous  to  represent  the  behaviour  of  both  of  the 
phases  in  single-phase  flow  by  a  single  line.  Lines  on  the 
same  plot  for  various  constant  gas-liquid  volume  ratios  can  then 
be  used  to  represent  the  two-phase  flow  behaviour.  The  shape 
of  these  lines  depends  on  the  flow  rates.  This  may  be  done 
by  defining  a  friction  factor  and  a  Reynolds  number  based  on 
the  mixture  properties.  Teletov  (15)  has  defined  a  friction 
factor,  f^pm*  and  a  Reynolds  number,  Rem,  in  this  manner. 

There  is  no  real  physical  basis  for  defining  the  friction 
factor  and  Reynolds  number  in  this  manner.  It  does,  however, 
allow  the  representation  of  the  single-phase  gas  and  single¬ 
phase  liquid  flow  by  a  single  line  on  R^  coordinates. 

Figures  13  and  14  are  such  representations  of  the 
pressure  drop  data  for  the  one-inch  and  the  two-inch  diameter 
pipes  respectively.  This  type  of  plot  shows  how  the  deviation 
from  the  single-phase  line  increases  as  the  liquid  rate  is 
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increased,  goes  through  a  maximum  and  then  decreases.  The 
deviation  from  the  single-phase  line  is  a  function  of  the 
Reynolds  number,  the  gas-liquid  volume  ratio  and  the  pipe 
diameter  for  the  fluids  investigated. 

Teletov  (15)  has  shown  this  deviation  to  be  a  function 
of  the  gas-liquid  volume  ratio  and  the  fluid  viscosity  ratio 
only.  Insufficient  data  are  presented  by  Teletov  to  allow 
a  critical  evaluation  of  his  work. 

Over  the  range  of  gas  rates  investigated  the  single¬ 
phase  friction  factor  is  fairly  constant,  especially  for  the 
one-inch  pipe.  This  makes  the  single-phase  pressure  drop 
directly  proportional  to  the  gas  density.  The  limited  two- 

i 

phase  data  taken  at  variable  density  clearly  indicates  that 
the  two-phase  pressure  drop  is  not  directly  proportional  to  the 
gas  density  or  the  mixture  density  at  approximately  constant 
superficial  liquid  and  gas  velocities.  It  is  shown  by 
Figure  15  that  the  effect  of  density  is  influenced  by  the 
liquid  and  the  gas  rates,  and  the  pipe  size.  In  the  annular- 
mist  flow  region  the  two-phase  pressure  drop  is  not  a  single 
valued  function  of  the  superficial  Reynolds  number  only.  This 

i 

makes  difficult  the  interpretation  of  the  data  of  other  inves¬ 
tigators  who  have  not  maintained  a  constant  mid-point  density. 
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The  experimental  work  for  this  investigation  was  conducted 
at  carefully  controlled,  constant  mid-point  pressure. 

The  manner  in  which  the  two  phases  orient  themselves 
in  the  pipe  has  a  great  effect  on  the  energy  loss  due  to 
irreversibilities.  The  higher  viscosity  fluid,  liquid  in 
this  case,  when  travelling  in  an  annular  film  against  the 
pipe  wal],  will  cause  a  greater  loss  of  energy  than  if  the 
same  amount  of  liquid  were  travelling  in  the  same  gas  stream 
as  a  dispersed  phase.  The  proportion  of  liquid  in  the  dis¬ 
persed  phase  increases  as  the  gas  velocity  is  increased.  This 
explains  the  greater  displacement  from  the  single-phase  pressure 
drop  of  the  two-phase  line  for  high  water  rates  at  low  gas 
rates.  The  effect  of  gas  velocity  on  the  film  thickness  seems 
to  be  more  pronounced  in  the  one-inch  pipe  than  in  the  two- 
inch  pipe.  Other  factors  which  would  be  expected  to  influence 
the  pressure  drop  are  the  shape  of  the  liquid  film,  distribu¬ 
tion  of  the  liquid  droplets  in  the  gas  phase,  the  impingement 
of  the  liquid  droplets  on  the  liquid  surface  and  formation 
and  collapse  of  the  surface  waves. 

The  comparison  of  these  data  with  those  of  other 
investigators  and  the  testing  of  various  correlations  was 
considered.  There  are  no  data  available  in  the  literature  for 
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the  range  of  variables  investigated  in  this  work,  making 
any  direct  comparison  impossible.  Considerably  more  work  is 
required  to  resolve  the  effect  of  the  physical  properties  of 
the  flowing  stream  and  make  possible  a  comparison  with  data 
obtained  for  different  systems.  The  correlation  of  Lockhart 
and  Martinelli  (7)  has  been  considered  by  other  investigators 
(7,  8)  and  has  been  found  not  suitable  in  this  region. 

The  analyses  of  Calvert  and  Williams  (17)  and  others 
(14,  15,  16,  18  and  19)  could  be  applied  to  these  data  to 
separate  the  water  flow  rates  in  the  film  and  the  gas  core. 

These  flow  rates  could  then  be  considered  separately  to  calculate 
the  energy  loss  due  to  irreversibilities.  This  work,  however, 
is  planned  as  a  separate  study. 

Further  investigations  are  required,  and  are  underway 
at  the  University  of  Alberta,  to  resolve  the  effect  of  gas 
density,  liquid  and  gas  properties  and  flow  rates  and  pipe 


size. 
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VII  CONCLUSIONS 

1.  A  mechanistic  approach  is  necessary  for  a 
better  understanding  of  the  two -phase  flow 
phenomenon . 

2.  The  pressure  drop  may  be  greatly  affected  by 
microscopic  changes  in  the  flow  pattern  within 
the  general  annular -mist  region. 

3.  Liquid  in  the  annular  film  causes  greater 

r 

energy  loss  than  the  same  liquid  in  the  gas  core. 

4.  Very  small  addition  of  water  to  a  high  velocity 
gas  stream  causes  a  large  increase  in  pressure 
drop  by  forming  a  liquid  film  on  the  inside 
pipe  wall. 

5.  At  higher  water  rates  the  proportion  of  the 

water  in  the  film  is  higher  for  lower  gas  velocities. 

6.  All  investigations  should  be  conducted  wherever 

> 

possible  at  constant  mid-point  conditions. 


NOMENCLATURE  * 


cross-sectional  area  of  the  pipe. 

sq.  ft. 

inside  diameter  of  pipe. 

ft. 

friction  factor , 

dimensionless 

dimensional  conversion  constant. 

lb. (mass)  ft. 
lb. (force)  sec2 

length  of  the  test  section. 

ft. 

pressure. 

lbs.  per  sq.  in 

Reynolds  number. 

dimensionless 

temperature  of  fluid  at  average 
conditions. 

°R 

pressure  drop  per  unit  length. 

lbs . 

ft?  perft- 

pressure  gradient  at  the  mid-point 

lb •  per  ft. 
ft? 

ratio  of  the  two-phase  pressure 

dimensionless 

gradient  to  the  single-phase  pressure 
gradient  if  the  same  quantity  of  gas 
was  flowing  alone  in  the  pipe  at  the 
same  mid-point  conditions. 


input  gas -liquid  mass  ratio. 

dimensionless 

input  gas -liquid  volume  ratio 
at  mid-point  conditions. 

dimensionless 

*  Special  terms  have  been  defined  in  their 
proper  context. 


61 


v 


V 


P 

M* 


specific  volume  of  the  liquid  at  cu.  ft.  per  lb 
mid-point  flow  conditions* 


superficial  velocity  of  fluid  at 
mid-point  conditions  based  on  total 
pipe  cross-section* 


ft.  per  sec. 


density  of  fluid  at  mid-point 
conditions 


lbs.  per  cu.  ft 


absolute  viscosity  of  the  fluid  lbs.  per  ft.  sec. 
at  average  mid-point  conditions 


Subscripts 


g  -  refers  to  the  gas  phase  or  values  calculated 
using  the  gas-phase  properties 

1  -  refers  to  the  liquid  phase  in  general  or  values 

calculated  using  the  liquid-phase  properties 

w  -  refers  to  the  water  or  values  calculated  using 
the  water  properties 

tp  -  refers  to  the  two-phase  flow 


tpg  -  refers  to  two-phase  flow  but  calculated  using 
gas  properties 

tpm  -  refers  to  two-phase  flow  but  calculated  using 
mixture  properties 

1*2, 3* 4, 5  -  refer  to  the  pressure  taps  on  each  test 

section  starting  at  the  upstream  end 


' 
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APPENDIX  A 


(1)  Composition  of  the  Gas 

(2)  Summary  of  Data  Describing  Test 
Section. 
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APPENDIX  A 

TABLE  A-l 

COMPOSITION  OF  THE  GAS 

Average  composition  of  the  gas  on  dry  basis  as  supplied 
by  the  Northwestern  Utilities 


Component 

Mole  Percent 

Methane 

86.18 

Ethane 

7.49 

Propane 

3.67 

i- Butane 

0.23 

n-Butane 

0.39 

i- Pentane 

0.17  i 

Hydrogen  Sulphide 

0 

Nitrogen 

1.58 

Carbon  Dioxide 

0.29 

Gravity  =  0.647  (Air  =  1.0) 

Average  gravity  was  recorded  for  each  day  and  it 


varied  between  0.64  and  0.66. 
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Water  Content 

Water  content  of  the  gas  was  measured  by  adsorbing 
the  water  from  a  measured  quantity  of  gas  on  anhydrous 
magnesium  perchlorate  contained  in  a  desiccant  tube. 
The  affinity  of  the  desiccant  for  hydrocarbons  is 
easily  satisfied  by  passing  a  small  amount  of  natural 
gas  (1  S.C.F.*)  through  it.  The  desiccant  is  capable 
of  absorbing  30  percent  of  its  weight  of  water.  Water 
content  of  the  gas  varied  between  0.000724  to  0.000806 
cm?  per  S.C.F. 


*  S.C.F.  means  a  cu.ftof  gas  at  14.4  psia  and  60°F. 
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TABLE  A- 2 


APPENDIX  A 


SUMMARY  OF  DATA  DESCRIBING  TEST  SECTIONS 


Measurement 


One-Inch  Pipe  Two-Inch  Pipe 


Pipe  Type 


Schedule  40  Schedule  40 

Aluminum  Pipe  Aluminum  Pipe 


Outside  Diameter, 

Inches  1.315  2.375 

Wall  Thickness, 

Inches  0.133  0.154 

Inside  Diameter, 

Inches  1.049  2.067 

Length  of  Upstream 

Calming  Zone,  ft.  24  22 

Length  of  the  Test 

Section,  ft.  40  80 

Number  of  Pressure 

Taps  5  i  5 

A  r- 

Distance  Between 

Pressure  Taps,  ft.  10  20 

Length  of  Downstream 
Calming  Section,  ft. 


15 


20 
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APPENDIX  B 


1 

2 

3 

4 

5 

6 

7 

8 

9 


One-Inch  Pipe  Experimental  Data 

Two -Inch  Pipe  Experimental  Data 

One-Inch  Pipe  Calculated  Data 

Two-Inch  Pipe  Calculated  Data 

One-Inch  Pipe  Pressure  Gradient  Ratio  Data 

Two-Inch  Pipe  Pressure  Gradient  Ratio  Data 

One-Inch  Pipe  Calculated  Data  Based  on 
Mixture  Properties 

Two-Inch  Pipe  Calculated  Data  Based  on 
Mixture  Properties 


Calculation  of  Data 
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SECTION  B-9 
CALCULATION  OF  DATA 


1.  Gas  Flow  Rate 

The  orifice  flow  equation  may  be  written  in  the 
form 

Qg  *  3600  *f  (17) 

where  Q  «  gas  flow  rate  in  standard  cubic  feet 

O 

per  second  (S.C.F.S.)  at  14.4  psia  and 
60°F. 


=  differential  pressure  in  inches  of  water. 
Pf  =  static  pressure  in  p.s.i.a. 

C'  =  F^  x  Fr  x  Y  x  Fpb  x  Ftb  x  Ftf  x  Fg  x 

FpV  x  Fm. 

The  values  of  all  the  factors  F^  ,  Fr,  etc.  are  obtained 
from  the  tables  in  Section  V  of  American  Gas  Associa¬ 
tion,  Gas  Measurement  Committee  Report  No.  3. 


2.  Gas  Density  at  the  Mid-point 

The  gas  law  may  be  written  in  the  form 


8 


29  G  P3 


Z  R  T 


(18) 
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where  P 

O 

G 

P3 

T 

R 

Z 


=  density  of  the  gas  at  the  mid-point 

o 

temperature  and  pressure,  lbs.  per  ft. 
=  gravity  of  the  gas  (Air  =  1.0) 

=  mid-point  pressure,  p.s.i.a. 

=  average  gas  temperature,  °R. 

•  gas  constant  =  10.73 
=  compressibility  factor 


3.  Superficial  Gas  Velocity  at  the  Mid-Point 

The  continuity  equation  may  be  written  in  the  form 


where 


fLSs 

ps  A 


(19) 


superficial  gas  velocity,  ft.  per  sec. 

gas  density  at  standard  conditions 

29  G  (14.4) 

"lx  10.73  x  520 


=  0.0488  lbs.  per  ft? 

Qg,  Pg  as  previously  defined 

o 

A  =  transverse  internal  pipe  area,  ftf 

Equation  (19)  reduces  to 

Vg  =  8.133  for  1-^  one-inch  pipe  (20) 

§ 

V„  =  2.098 

S  pg 


for  the  two-inch  pipe  (21) 


- 
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4.  Water  Rate 

The  orifice  flow  equation  for  incompressible  fluids 
may  be  written  in  the  form 

c"  y nr  (22) 

orifice  constant 

differential  pressure 

2 

(Instrument  Reading)  x  8 

The  flow  rates  were  directly  obtained  from  the  orifice 
calibration  curves.  Figure  C-l. 


Qw  - 
c"  - 

K  - 


5.  Superficial  Water  Velocity 


V  ■  Qw 
w  A 


(23) 


where  Vw  -  superficial  water  velocity,  ft.  per  sec. 


6.  Reynolds  number  Based  on  Mixture  Properties 

Pt 


Re, 


m 


■  •[ 


pg  Vg 

*g 


w  Vw~l 

Hw  J 


(14) 


where  Rem 


Reynolds  number  based  on  mixture 
properties . 

gas  viscosity,  -57^ — 

0  ft.  sec. 

0.679  x  10*5 

lb. 

water  viscosity, 

0.672  x  10~3 

gas  density  at  the  mid-point,  lbs.  per 

ft?  -  0.47. 


Ill 


Pw  =  water  density  at  the  mid-point,  lbs. 
per  ft?  =  62.43. 

D  =  internal  pipe  diameter,  ft.*. 0874  and 
.1722. 

V  ,  V  as  previously  defined 
g  w 


7.  Friction  Factor  based  on  input  Mixture  Properties 


'tpm 


/dP\ 

^dlP  gc  D 

2  (Vg  +  Vw) 


(13) 


m 


where 


tpm 


friction  factor  based  on  mixture 


properties. 


(<£) 

ML' 


pressure  gradient  at  the  mid-point. 


*  32.1740 


^m 

lbf 


ft. 

sec? 


m 


=  mixture  density 

=  pg  +  pw  Vw 

V  +  V 
vg  +  vw 


D,  Vg,  Vw,-  Pg,  Pw  as  previously  defined. 


8.  Reynolds  Number  based  on  Gas  Properties 


where 


Res 

Reg 


D  pg  vg 

Reynolds  number  based  on  gas  properties 
at  the  mid-point. 


- 
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The  equation  reduces  to 


Reg  - 

12878 

pg  vs 

for 

the 

one -inch 

pipe 

(24) 

Reg  = 

25377 

pg  vg 

for 

the 

two -inch 

pipe 

(25) 

9.  Friction  Factor  based  on  Gas  Properties 


/dpx 

£  —  (<«  )  §c  ® 

ttpg  “  -a—~ - - - 


2  V 


(26) 


g 


where 


tpg 


friction  factor  based  on  gas  pro¬ 
perties. 


The  equation  reduces  to 


'tpg 


^  €>  z  7? 

Kg  vg 

for  the  one-inch  pipe 


(27) 


tpg 


2.777  (£) 

dL  pg  V 


g 


(28) 


for  the  two-inch  pipe 
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APPENDIX  C 

CALIBRATION  OF  THE  EQUIPMENT 

C-l  Gas  Orifice  Meter 

C-2  Pneumatic  Recorders 

C-3  Pneumatic  D/P  Cells 

C-4  Orifice  Calibration  for  the  D/P  Cell 

C-5  Electrical  D/P  Cell  and  Recorder 
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APPENDIX  C 

CALIBRATION  OF  THE  EQUIPMENT 

C-l 

Gas  Orifice  Meter; 

The  static  pressure  element  on  the  gas  flow  recorder 
was  calibrated  against  a  dead-weight  pressure  gauge.  No 
correction  to  the  chart  reading  was  necessary  over  the  range 
of  interest.  The  differential  element  was  calibrated  against 
a  mercury  manometer  and  no  correction  was  necessary  to  the 
chart  reading. 

Individual  orifices  were  not  calibrated.  The  calibra¬ 
tion  data  of  American  Gas  Association*  Gas  Measurement  Committee 
Report  No.  3*  were  used  for  all  flow  calculations.  A  few  spot 
checks  were*  however*  made  on  two  of  the  orifices  by  the  use 
of  a  precalibrated  critical  flow  prover.  In  all  cases  the 
gas  measurement  by  the  orifice  meter  was  within  1%  of  the 
measurement  by  the  critical  flow  prover. 


C-2 

Pneumatic  Recorders ; 


Pneumatic  recorders  were  used  to  record  the  pressure 
drop  in  the  integral  orifice  D/P  cell  for  water  measurement. 
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mid-point  pressures  in  the  two-inch  and  the  one-inch  test 
sections,  and  the  pressure  drop  in  the  two-inch  test  section. 
The  recorders  were  calibrated  against  a  test  gauge  by  supplying 
a  pressure  of  3  to  15  psig  to  the  recorder. 


C-3 

Pneumatic  D/P  Cells ; 

The  pressure  across  the  diaphragm  of  a  D/P  cell  was 
measured  on  a  manometer  and  the  output  of  the  D/P  cell  was 
recorded  on  a  precalibrated  recorder.  The  range  of  all 
pneumatic  D/P  cells  was  checked  in  this  manner,  including 
the  integral  orifice  D/P  cell  for  which  the  orifice  was  re¬ 
moved  for  range  calibration. 

C-4 

Orifice  Calibration  for  the  D/P  Cells 

Orifices  for  the  integral  orifice  D/P  cells  were 
calibrated  by  collecting  the  water  for  certain  periods  of 
time  and  weighing  it.  The  calibration  data  for  0.25,  0.0595, 
0.034,  0.0995  and  0.02  inch  orifices  is  plotted  in  Figure  C-l 
as  instrument  (recorder)  reading  versus  water  rate  in  cm?/sec. 
The  recorder  had  a  square  root  chart  and  a  range  of  0-800 


inches  of  water. 
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05 

Electrical  D/P  Cell  and  Recorder; 

The  electrical  recorder  and  the  electrical  D/P  cell 
were  checked  in  the  same  manner  as  the  pneumatic  instruments. 
They  were  both  found  to  be  more  accurate  than  the  pneumatic 
instruments . 


WATER  RATE  ,  CITT  per  See. 


water  rate  ,  cm.  per  sec. 

FIGURE  C-l  ,  CALIBRATION  OF  WATER  D/p  CELL  ORIFICES 
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